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Abstract 

Five-dimensional cosmological models with two 3-branes and with a buck cosmological constant 
are studied. It is found that for all the three cases (A = 0, A > 0, and A < 0), the conventional 
space-time singularity "big bang" could be replaced by a matter singularity "big bounce" , at which 
the "size" of the universe and the energy density are finite while the pressure diverges, and across 
which the universe evolves from a pre-existing contracting phase to the present expanding phase. 
It is also found that for the A > case the brane solutions could give an oscillating universe model 
in which the universe oscillates with each cosmic cycle begins from a "big bounce" and ends to a 
"big crunch", with a distinctive characteristic that in each subsequent cycle the universe expands 
to a larger size and then contracts to a smaller (but non-zero) size. By studying the gravitational 
force acted on a test particle in the bulk, a gravitational stability condition is derived and then is 
used to analyze those brane models. It predicts that if dark energy takes over ordinary matter, 
particles on the brane may become unstable in the sense that they may escape from our 4D-world 
and dissolve in the bulk due to the repulsive force of dark energy. 
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I. INTRODUCTION 



In the brane-world scenarios, our conventional universe is a 3-brane embedded in a higher 
dimensional space. While gravity can freely propagate in all dimensions, the standard matter 
particles and forces are confined to the 3-brane only [1]. In recent years, five- dimensional 
(5D) brane cosmological models have received extensive studies [2-9]. It is noticed that one 
of the many interesting features of brane models concerns the big bang singularity: Due to 
the existence of extra dimensions, the conventional big bang singularity could be removed. 
So the big bang is perhaps not the beginning of time but a transition from a pre-existing 
phase of the universe to the present expanding phase, and our universe may have existed for 
an infinite time prior to the putative big bang. In the ekpyrotic model [10], it was suggested 
that the universe was produced from a collision between our brane and a bulk brane. In the 
cyclic model [11], the universe undergoes an endless sequence of cosmic epochs that begin 
with a big bang and end in a big crunch. In the bounce models [7-9,12-14], it was shown 
that the scale factor could evolve across a finite (but non-zero) minimum which represents a 
bounce (as opposed to a bang). The purpose of this paper is to study the bounce property 
of brane models more generally. 

In a previous paper [15], a five- dimensional big bounce cosmological solution with a non- 
compact fifth dimension was presented. By using this solution as valid in the bulk, a global 
brane model is derived in Ref. [9], in which the model has two 3-branes with the extra 
dimension compactified on an Si/Z 2 orbifold. This brane model is of the type of Binetruy, 
Deffayet and Langlois [2,3] in which no cosmological constant was introduced in the bulk. 
In this paper, we are going to generalize it by adding a cosmological constant in the bulk, 
and then to study evolutions of the models. 

The plan of this paper is as follows. In Section II, we look for general solutions for 
cosmological models with two 3-branes and for three cases with A = 0, A > and A < 0, 
respectively. In Section III, we study the gravitational force acted on a test particle in the 
vicinity of a brane and derive a stability condition. In Section IV, we give several simple 
exact solutions as an illustration, and study the global evolutions and the stability of the 
brane models as well as the big bounce singularity. Section V is a short discussion. 
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II. GENERAL BRANE SOLUTIONS WITH A COSMOLOGICAL CONSTANT 



Let the 5D metric being 

dS 2 = B 2 dt 2 - A 2 



dr 2 



+ r 2 dn 2 



dy 2 



(1) 



1 — kr 2 

where B = B(t y) and A = A(t t y) are two scale factors, k(= ±1 or 0) is the 3D curvature 
index, and dVL 2 = (d9 2 + sin 2 Odcp 2 ). This metric describes 5D cosmological models with 
spherical symmetry in 3D and with a static fifth dimension. Now we consider brane models 
with the extra fifth dimension y compactified on a small circle for which the first brane is 
at y = y\ = and the second is at y = y% with y% > 0. The 5D Einstein equations read 

Gab = k 2 (Ag AB + T AB ) 

Tab = [0°* + P*) ~ Pi9 aP ] 9aA9/3B$ (y ~ Vi) , (2) 
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where upper case Latin letters denote 5D indices (0,1,2,3;5), lower case Greek letters denote 
4D indices (0,1,2,3), k 2 is the 5D gravitational constant, u a = dx a /ds is the 4D velocity in 
comoving coordinates with ds being the 4D line-element and u a = (u°, 0,0,0), and and 
Pi are energy density and pressure on the z-th brane, respectively. With use of the metric 
the non- vanishing equations in (j2J are found to be 
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where an overdot and a prime denote partial derivatives with respect to t and y, respectively. 

Note that there are different ways to express the solutions of above equations (see, for 
example, [3]). Here we follow our previous work [9,15,16] for convenience. We see that 
equation (J3J) can be integrated once, giving 

where is an arbitrary function. Using this relation to eliminate B from equation (JHJ), 
we obtain 
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A 2 = 2(^ 2 + k) . (8) 



Also, with use of the relation (J7J), equation (JHJ) can then be integrated, giving 

+ k _ a' 2 ) A 2 = -k 2 AA a + F , (9) 
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where F is an integration "constant" which actually could be a function of y. 

In the bulk, the two equations (JHJ) and (JHJ) require that F = K = const. So the Einstein 
equation in the bulk is 

(A 2 )' 2 + ^k 2 AA 4 = AA 2 [^ 2 {t) +k] -4=K . (10) 

Meanwhile, we can verify that equation (jUJ) is then satisfied identically. 

From a known exact bulk solution of the equation (jlOj) we can easily extend it from bulk to 
branes and obtain a global two-brane solution. The S\/Z 2 symmetry requires that we firstly 
should write a bulk solution from gAB(t,y) to gABit, \y\). Then, according to Israel's jump 
conditions, the two scale factors A(t, \y\) and B(t, \y\) are required to be continuous across 
the branes localized in y — y{. Their first derivatives with respect to y can be discontinuous 
across the branes. And their second derivatives with respect to y can give a Dirac delta 
function. Thus the resulting terms with a delta function appearing in the LHS of equations 
(JHJ) and (JHJ) must be matched with the corresponding terms containing a delta function in 
the RHS of them in order to satisfy the field equations. For we have, for instance, 

A , :^A^+[A'] i 8(y-y i ), (11) 

where A" represents terms in A" that are not contributed to the delta function 8(y — y-i), 
and [A'] i is the jump in the first derivative across y = y^ defined by [A'] i = A\yf) —A'(y^). 



The Z2 reflection symmetry leads to 
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Thus the Einstein equations (J3j) and © on the branes give 
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Meanwhile, the 5D conservation law Tj±. B = gives the usual 4D equation of conservation 
on the branes: 



Pi + 3 (pi+Pi) ^ = 



(14) 



We also need to define the Hubble and deceleration parameters on the branes appropri- 
ately. Be aware that the proper time on a given y =constant hypersurface is defined by 
dr = Bdt, so, with use of the relation (jZJ), the Hubble and deceleration parameters are 
defined as (see [9]) 
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So, on the branes, we have 
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Using ([TBI) and ([T3l in (fTOjl . we obtain 
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Thus we recover the induced Friedmann equation on the branes that was discussed and 
studied widely in literature [2-6,9,17]. In what follows we consider the three cases for the 
brane models with A = 0, A > 0, and A < 0, respectively. 



A. TYPE I BRANE MODELS: A = 



For A = 0, the general solutions of (fTUjl are found to be [9] 



A 2 =^ + k)y 2 -2u\y\ + ^. (18) 
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This type of solutions contains two arbitrary functions n(t) and v(t) and two constants k 
and K. Differentiating this A 2 with respect to y, we obtain 

d\y\ 
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dy 



(19) 



Therefore we get 
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where the first brane is at y = y\ = and the second is at y = y 2 > 0. Then, from (J7j) we 
get 
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So equation (fTHj) gives 
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B. TYPE II BRANE MODELS: A > 

For A > 0, the general solutions of (fTUj) can easily be found. Then, by changing y to \y\, 
we obtain 



A 2 = a(t) cos (9(t) + L~ x |y|) + 2L 2 (/i 2 + k) , 

L- x = 



a(t) = 2Ly L 2 (fx 2 + k) — K 



2k 2 A 



(24) 



6 



where //(£) and 8(t) are two arbitrary functions. By differentiating A 2 with respect to y, we 
obtain 

r\ I I 

2 A A' = -L- X a{t) sin (0 + L' 1 \y\) -jp- . (25) 



So we have 
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Using these and the relation B = Aj /x, we get 
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So equation (JT3)l gives 
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C. TYPE III BRANE MODELS: A < 

For A < 0, the general solutions of (fTUj) are found to be 

A 2 = a(t) cosh (f (t) + L- 1 \y\) - 2L 2 (u 2 + k) 



a(t) = 2L^L 2 {tf + kf + K , L~ l = J ~k 2 A , (30) 

where k and K and L are constants, /x = /x(i) and £ = £(£) are two arbitrary functions. 
Differentiating A 2 with respect to y, we obtain 
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So we find 
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Then, using these and the relation B — A/ ft, we find 
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So the equation (fTHj) gives 
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III. A GRAVITATIONAL STABILITY CONDITION FOR PARTICLES ON 
BRANES 



Generally speaking, brane models require the fifth dimension to be compactified to a small 
size 2/2. So a stable brane model should have a stable size. The well-known Goldberger-Wise 
mechanism [22] was to use a scalar field in the bulk to stabilize the size of the fifth dimension. 
In this way, the model has the tendency to recover it's size after a perturbation. Here, a 
perturbation means a small change for the size. 

Here we wish to say that even if the size of the fifth dimension is somehow stabilized, one 
still can ask question whether particles on the brane may leave the brane and escape into 
the bulk. Arkani-Hamed et al [1] pointed out that in sufficiently hard collisions, particles on 
the brane can acquire momentum in the extra dimensions and escape from our 4D world, 
carrying away energy. If this happens continuously, it will cause another kind of instability 
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problem. Note that here the instability just means brane particles are not in a stable position 
along the fifth dimension. Now a natural question is that once a particle left the brane, will 
it escape in the bulk forever or return to our brane again? 

To answer this question we noticed that the bulk discussed in this paper only contains 
a cosmological constant term A and so is empty. Therefore, it is reasonable to expect that 
particles inside the bulk should obey the 5D geodesic equation 

which was used previously [18] for a similar purpose as here. From this equation, a 5D 
gravitational force acted on the bulk test particle can be defined as 

a „ a doc doc 



BC dS dS ' 



(37) 



For simplicity, we assume that the particle is temporary at rest in the bulk. Then, along 
the fifth dimension, (jS7j) gives 

F 5 = -f • (38) 

In the vicinity of the i-th brane, we use the general results (fT3*|) and (|12jl in ()38j1 . Then we 
obtain 

F\vt) = -F\v7) = -y + 5ft) • (39) 

If this force F 5 is attractive, i.e., F b {yf) < and F 5 (y^) > 0, the particle would be "pulled" 
back to our brane. Thus we obtain a new kind of stability condition for the i-th brane as 

Pi + \vi > . (40) 

This is a reasonable condition which holds for ordinary matters including dark matter. 
However, if the brane matter contains a dark energy component such as a cosmological 
constant or quintessence, the condition (|4T)j) may or may not hold, depending on how much 
the dark energy is contained on the brane. For example, we let 

pi = pi + A; , Pi=pi-\i , (41) 

where \ is a cosmological term on the i-th brane. Then the condition (|4T)|) becomes 

3 1 - 

Pi + 2& > 2 l% ■ (42) 
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Suppose initially condition (}4*2"j) holds, then, as the universe expands, both p, and pi decrease 
while Aj keeps unchanged. So gradually the universe will enter in an unstable stage in which 
particles and energy on the brane may escape in the bulk and the 4D conservation law of 
energy (|14j) may not hold anymore. Recent observations [19,20] reveal that presently we are 
living in an accelerating stage dominated by a dark energy term such as \ in (|42j) . Thus we 
see that both the acceleration of our universe and the unstable nature of the brane particles 
could be explained as due to the same repulsive force of dark energy. In the following section 
we will show that some brane models do not satisfy the condition (J4~U|) . 

IV. SOME SIMPLE BIG BOUNCE BRANE MODELS 

In Section II we have obtained three types of global brane solutions corresponding to 
A = 0, A > and A < 0, respectively. Each type contains two arbitrary functions: p(t) 
and v(t) for type I, p(t) and 9(t) for type II, and pit) and for type III. All the two 
scale factors A(t,y) and B(t,y) and the densities pi and pi on the branes are expressed via 
the two functions. From the relation B = A/p and metric we see that the form of Bdt 
is invariant under an arbitrary coordinate transformation t — > t(t). This freedom could be 
used to fix one of the two arbitrary functions. Another freedom may correspond, as is in 
the standard general relativity, to the unspecified equation of state of matter. So, generally 
speaking, if the matter content on the first brane is known, then the two arbitrary functions 
could be fixed. Then the whole solutions could be fixed too. Then we will know the matter 
content on the second brane. Therefore, the brane solutions obtained in Sec. II are quite 
general. 

To compare these cosmological solutions with observations, we need know clearly the 
matter content on our brane. This may need careful analysis and might be complicated 
mathematically, and we are not going to do it here in this paper. However, as an illustration, 
we will pick up several simple exact models in the following and to exhibit typical features 
of brane cosmology. 
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A. A SIMPLE TWO-BRANE MODEL WITH A = 



Consider the A = brane solutions (|18p. for which we choose 

k = , K = 1, v(t) = t b /t , n(t) = (2ty 1/2 
where t b > is a constant. In this way, the solution (fTHj) becomes 
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On the second brane we have 
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4t 2 - (j/ 2 - 2t b y 4ti + (y 2 -2t b y 
It is easy to see for the first brane that there is a critical time t = t b > at which the scale 
factor A\(t) reaches to a non-zero minimum 2^/T b . After this critical time, A\(t) increases, 
implying the universe is expanding. For t t b , we have Ai(t) — > y/2~i, B\(t) — > 1, 
Pi — > 3pi, and the universe evolves as if in the radiation-dominated standard Friedmann 
model. Before this critical time, the universe contracts from infinity. It was shown [9,21] 
that as the coordinate time t varies from zero to t b , the proper time r varies from — oo to 
r b , where r b is a finite constant and corresponds to t b . Thus we can call t b as to represent 
a big bounce (as opposed to a big bang), and before this bounce the universe has existed 
for an infinitely long time. Meanwhile, rather than the big bang singularity, which should 
correspond to Ai(t) = 0, this big bounce singularity corresponds to Bi(t) = 0. From the 
solution (|45|) we can also see that as t varies across t b , the energy density p\ remains finite 
while the pressure p± changes from — oo to +oo . This implies a phase transition happened 
at the bounce. Here we plot the evolution of the scale factor Ai(t) with t b — 1 in Figure 1. 
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FIG. 1: Evolution of the scale factor Ai(t) = ^2(t + l/t) in a A = brane cosmological model. 



The evolution of the second brane is similar as the first one with the bounce point at 
t = U2 = \h — 2/2/2 1 . From equation (jlBj) we see that whether the energy density on the 
second brane is positive or negative depends on the "size" of the fifth dimension y 2 . If 
j/2 > 2tfe, we have p 2 > 0. If y 2 < 2tb, we have p 2 < 0. 

Now let us look at the stability condition pHjl . Using PHJ) and (jjHj) in we obtain 
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For the first brane we find that pi + |pi > for t > t 6 . So, in the expanding stage (after the 
bounce), the first brane is stable. There are two cases for the second brane. If y 2 > 2tb, we 
have p 2 > and p 2 + jp 2 > in the stage t > U2 where U2 = 1 2/2 — 2tb\ represents the bounce 
time. So after the bounce this brane is stable. If y 2 < 2%, we have p 2 < and p 2 + |p 2 < 
in the stage t > t^- So after the bounce this brane is unstable. Thus, to obtain a model 
with two stable branes, we must have y 2 > 2%- An interesting special case is y 2 = 4tb for 
which we have p 2 = pi > and the model is completely symmetric and particles on the 
branes are stable. 
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B. AN OSCILLATING UNIVERSE MODEL WITH A > 



The notion of an oscillatory universe can be traced back to 1930's and has continued to 
attract interest [23-25,11]. Tolman [23] discussed it within the framework of general relativity 
assuming a closed universe (k = +1) in which the universe undergoes a sequence of cycles of 
expansion and contraction. Dicke and Peebles [24] restudied Tolman's model and pointed out 
that an oscillating universe could provide an escape from some of the cosmological problems 
such as the horizon and homogeneity problems. u As the universe ages" , they wrote, "more 
and more of it becomes visible that earlier was beyond the horizon and presumably causally 
disconnected from us .... How then are we to understand the remarkable familiarity of 
the objects just appearing on the horizon? Perhaps by tracing the evolution back through 
the big bang to an earlier collapsing phase" . However, Tolman's oscillatory models were 
constrained by having to pass through the big bang singularity in which the energy density 
and temperature diverge. Therefore, Dicke and Peebles expected that "some future and 
better theory might show that the collapse of the universe would lead to a 'bounce ' instead of 
a singularity" . 

Now let us consider the type II (A > 0) brane solutions. As we discussed at the begin- 
ning of this section that the metric form of Bdt is invariant under an arbitrary coordinate 
transformation t — > t(t). Meanwhile, the scale factor A(t,y) in (|2*%Jl contains two arbitrary 
functions fi(t) and 9(t). So, by choosing the time coordinate properly, we can set 6{t) = —ixt 
without loss of generality. So the cosine term in the solution leads generally to an oscillating 
universe model. For illustration, we choose 



fi 2 (t) = t, 



0(t) = — Tit , 
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So on the branes we have 
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FIG. 2: Evolution of the scale factor Ai(t) = \J2t + 2Vt 2 - 1 cos (vri) in a A > oscillatory brane 
cosmological model 

We plot the evolution of the scale factor A\(t) in Figure 2. From this figure we see clearly 
that the big bang spacetime singularity of the standard cosmology is replaced here by a 
series of smooth big bounces at t — i& ~ 1, 3, 5, 7, .... 

Another typical feature of Figure 2 is that there is a "beginning of time" at t = 1 
in the model and the oscillating amplitudes of the followed cycles increase monotonously. 
This reminds us of the well-known entropy problem of the old oscillatory universe model. 
Tolman [23] pointed out that if the total entropy of the universe can only increase, then, 
in an oscillatory model, the entropy produced during one cycle would add to the entropy 
produced in the next cycle, causing the oscillating amplitude of each cycle to be larger than 
the one before it. Extrapolating backward in time, the universe would have to be finite 
cycles old. More discussions about the entropy problem can be found in literature [24,11]. 
We see that our oscillating universe model described by equation (|5T|) and Figure 2 coincides 
with Tolman's description perfectly and, therefore, could provide a realized framework to 
discuss the entropy problem. 

The two equations (f2"K|) and (|2"Hj) become 
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So the energy density on the branes changes sign periodically with time, implying that 
the universe may have a negative energy density. This is an unusual feature of the oscillatory 
brane model. It is reasonable to assume that the size the fifth dimension is much smaller 
than the period of each cycle of the universe, i.e., y 2 <C 2tt. Then from equations (JoTj) - (J53J) 
we find A 2 (t) « Ai(t), p 2 ~ —pi and p 2 ~ —pi- So if our brane has a positive energy density 
at present stage of the universe, the other brane would have a negative energy density. 
Applying this in the stability condition (j4*Uj) . we find that (p 2 + |p 2 ) ~ — (pi + |Pi)- This 
means that if one of the two branes is at a stable stage, the other one must be at an unstable 
stage. As a whole, we conclude that brane particles in this oscillating model are not stable. 



C. A SIMPLE UNIVERSE MODEL WITH A < 

For this type of brane solutions ()30p. we take a similar choice as in equations 

p 2 (t) = t, £(t) = -t , 

L = K = 1, k = . (54) 

Then we obtain 



A 2 = 2 Vt 2 + 1 cosh {t - \y\)-2t 



A\ = 2Vt 2 + 1 cosh (t - - 2t . (55) 

We plot the evolution of Ai(t) of the first brane in Figure 3, from which we see that Ai(t) 
reaches to its minimum at the bounce point t = tj, « 0.5. Before the bounce, the universe 
contracts from infinity; after the bounce, the universe expands to infinity again. 
From equations ()34|) and (|35|) we obtain 
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" R = "^«l — m — a\ — • (56) 
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FIG. 3: Evolution of the scale factor A\(t) = \J lyj t 2 +T cosh t — 2t in a A < brane cosmological 
model. 

and 



.2 



« 92 



6 ^t 2 + 1 sinh (t-2/ 2 ) 



2 2 9 /-y/F+lsinh^-^A 4 v / ^TTsinh(t-y 2 ) 

KP2 = p ^ * J + a! ■ (57) 

Thus the energy density is positive on the first brane and negative on the second brane. If 
we assume \t\ 3> 2/2, then we have P2 ~ —p\ and P2 ~ — Pi- So probably we are living on the 
first brane. As for the stability condition (jH3j) , we also have (p 2 + ~ — (Pi + |pi) as 
in the above oscillating model. So particles in this model are also unstable. 



V. DISCUSSION 



In this paper we have derived, in Sec. II, a class of five-dimensional cosmological solutions 
with two 3-branes and with the fifth dimension being static and compactified on a small 
circle. The bulk contains only a cosmological constant A. It is found that for all the three 
cases of A (A = 0, A > 0, A < 0) the solutions contain two arbitrary functions of time. 
One of these two freedoms might be explained as due to the unspecified time coordinate in 
the 5D metric leaving another to account for various contents of the cosmic matter. In 
Section III we have used the 5D geodesic equations to study the gravitational force acted 
on a test particle in the vicinity of a brane. This force could be interpreted as generated 
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by matters on the brane. By requiring this force being attractive and so to grip particles 
from escaping into the bulk, we have derived a physically reasonable gravitational stability 
condition as given in equation (|4U|) . For illustration and for simplicity we presented three 
simple exact models in Section IV by choosing the two arbitrary functions properly. From 
these simple models we found that the conventional space-time singularity "big bang" could 
be replaced in brane models by a matter singularity "big bounce" at which the "size" of the 
3D space is finite and the energy density does not diverge, while the pressure diverges. This 
enable us to expect that in brane cosmological models the "history" of our universe could 
be traced back across the big bounce to a pre-existing phase. This is clearly of great interest 
and deserve more studies. The stability of brane particles of these simple models are also 
studied. 

Here we want to discuss more about the oscillating universe models given in equations 
()50|) - ()53|). As pointed out by Tolman that an oscillatory model could resolve the horizon 
and the homogeneity problems. However, the main difficulty of Tolman's oscillatory model 
is having to pass through the big bang space-time singularity. Now brane models could 
remove the big bang singularity in a satisfactory way and thus rescued Tolman's old model. 
Meanwhile, Tolman's entropy problem also get resolved. 

We should emphasis that the present work is exploratory. Be aware that the general 
brane solutions given in this paper contain two arbitrary functions. This would enable us 
to discuss more observations such as the acceleration of the universe [23]. We leave these 
studies in the future. 
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